Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 16 October 2008 



(MN m^K style file v2.2) 



Implications of Recent Measurements of the Milky Way 
Rotation for the Orbit of the Large Magellanic Cloud 



oo 
O 

o 

(N 
-t— » 
o 

O 

in 



43 

6 



Genevieve Shattow and Abraham Loeb 

Harvard-Smithsonian Center for Astrophysics, 60 Garden St, Cambridge, MA 02138 



16 October 2008 



ABSTRACT 

We examine the implications of recent measurements of the Milky Way (MW) rotation 
for the trajectory of the Large Magellanic Cloud (LMC). The ~ 14 ± 6% increase 
in the MW circular velocity relative to the IAU standard of 220 km s" 1 changes the 
qualitative nature of the inferred LMC orbit. Instead of the LMC being gravitationally 
unbound, as has been suggested based on a recent measurement of its proper motion, 
we find that the past orbit of the LMC is naturally confined within the virial boundary 
of the MW. The orbit is not as tightly bound as in models derived before the LMC 
proper motion was measured. 
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1 INTRODUCTION 

Recently, Kallivayalil et al. (2006; hereafter K06) mea- 
sured the proper motion of the the Large Magellanic Cloud 
(LMC), and pioneered a first assessment of its 3D veloc- 
ity vector. Based on this measurement, Besla et al. (2007; 
hereafter B07) concluded that the LMC was unlikely to have 
passed near the Milky Way (MW) before and was most likely 
formed outside of the boundaries of the Galaxy, in contrast 
to traditional scenarios (see references in van der Marel et 
al. 2002; hereafter vdM02). This new conclusion is intriguing 
given that there are no other examples of massive gas-rich 
galaxies like the LMC within the much bigger volume that 
separates the MW from its neighboring galaxy, Andromeda 
(M31). Other recent studies have also found the B07 results 
difficult to accept and propose alternate strategies of bind- 
ing the LMC and MW either through MOND gravity (Wu et 
al. 2008) or by giving the LMC and its smaller counterpart, 
the Small Magellanic Cloud, a common halo (Bekki 2008). 

A couple of years after K06 published their findings, 
Piatek et al. (2008; hereafter P08) confirmed independently 
their results to within one standard deviation, although at 
the lower end of the inferred range of values. Also, within 
the past five years, the circular velocity of the MW and 
the distance between the Sun and the galactic center have 
been updated by Reid & Brunthaler (2004; hereafter RB04) 
and Gillessen et al. (2008; hereafter GGE08), respectively. 
These increased the likely circular velocity of the MW from 
the IAU standard of V clrc = 220 km s" 1 to 251±15 km s~\ 
A value of 220 km s _1 now corresponds to a reduction of the 
best-fit value by two standard deviations (and equivalent to 
moving the Sun a total of 1 kpc closer to the center of the 
Galaxy). Uemura et al. (2000) used parallax measurements 



from Hipparcos and SKYMAP to obtain a similar value of 
V c i rc = 255±8 km s _1 . 

The rotation speed of the MW affects the analysis of 
the past LMC orbit in two ways. First, because the proper 
motion of the LMC is measured relative to the solar sys- 
tem which orbits the Galaxy, it is necessary to know the 
rotational velocity of the Sun in order to transform to the 
Galactocentric frame (see, e.g. vdM02). Second, the depth 
of the gravitational potential well of the MW (involving its 
estimated mass and scale radius) depends on the normaliza- 
tion of its rotation curve. B07 adopted the IAU standard in 
their analysis instead of the modified values for the Milky 
Way's rotation. In this Letter, we examine the implications 
of the change in the MW parameters for the LMC orbit (also 
with the updated P08 value). In the particular geometry of 
the LMC orbit, both of the above-mentioned effects make 
the LMC more gravitationally bound to the MW owing to 
an increase in Vdrc- Despite the small fractional magnitude 
of the correction in circular velocity (~ 14±6%) , we find that 
the qualitative nature of the LMC orbit changes. Instead of 
the LMC being possibly unbound as suggested by B07, the 
~ 10% decrease in the LMC velocity and the ~ 50% increase 
in the MW mass lead us to conclude that the LMC's past 
trajectory was probably confined within the virial radius of 
the MW. The apogalacticon distance of the orbit is com- 
parable to the MW virial radius, as expected for a satellite 
that had formed at the outer edge of the Galactic halo. 

Traditional studies of the LMC's orbit around the MW 
(e.g., Murai & Fujimoto 1980; Lin & Lynden-Bell 1982; Gar- 
diner et al. 1994; Lin et al. 1995; Gardiner & Noguchi 1996; 
vdM02; Bekki & Chiba 2005; Mastropietro et al. 2005; Con- 
nors et al. 2006) considered the MW as an isolated galaxy. 
While this might have been acceptable for studies where 
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the LMC's orbit was confined well within the MW halo, the 
high LMC velocity measured by K06 and P08 implies (B07) 
that the apogalacticon could extend beyond the edge of the 
MW's halo - where the gravitational influence of M31 is 
non-negligible (see Table 1 in B07). Thus, we also include 
the tidal effect of M31 in our calculations. 

In 32] we describe our adopted model for the mass dis- 
tribution of the MW halo. We then calculate the past LMC 
orbit (33| and the effect of M31 on it Finally, we discuss 
the implications of our results in SJU 



2 METHOD 

Following B07, we adopt a Navarro, Frenk, & White (1996; 
hereafter NFW) mass profile for the dark matter distribu- 
tion in the MW halo, and include dynamical friction (Chan- 
drasekhar 1943, Hashimoto 2003) in calculating the LMC 
motion through the MW halo. We also add the gravita- 
tional potential of M31 in the form of another NFW profile 
at a present-day Galactocentric distance of 780 kpc (Mc- 
Connachie et al. 2005; Cox & Loeb 2008 and references 
therein). The full gravitational potential <&tot as a function 
of radius r in each galaxy (either MW or M31) includes con- 
tributions from a disk (&disk), a bulge (<&buige), and a dark 
matter halo (&nfw), 



$tot(r) = <t>disk(r) + $bulge(r) + $NFw(r), 

where (Xue et al. 2008) 
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We adopt the standard ACDM cosmological parameters 
O m = 0.3, Ho = 70 km s^Mpc" 1 (Komatsu et al. 2008), 
and an overdensity of 8th = 340 (Wechsler et al 2002; Klypin, 
Zhao, & Somerville 2002), with 



47T 3 3 

M v i r = —^-pcrQm$thTvir QC ^eircj 



(6) 



and p C r = 3-ffo /8ttG. 

Some of our initial conditions (e.g., M vir and r„j r ) de- 
pend on the distance of the Sun from the center of the MW 
galaxy, as will be discussed in |3] but we consistently adopt 
M di3k = 4 x 1O 1O M , M bulge = 0.8 x 1O 1O M , a disk scale 
length b = 3.5 kpc, a disk scale height co = 0.7 kpc, and 
a halo concentration c = 12 for the MW; and Mdisk = 
7 x 10 10 M Q , M bu i a e = 1.9 x lO lo A/ , M vir = 1.6 x 10 12 M Q , 
b — 5.7 kpc, Co = 1.14 kpc, r„j r = 300 kpc, and c = 12 for 
M31, as suggested by Klypin et al. (2002). Aside from the 
revised normalization of M v i r and r v i r based on the mod- 
ified value of V c i r c, we have used the same mass profile as 
B07. 



3 CONSTRAINTS ON THE ORBITAL 
HISTORY 

Previous calculations of the history of the LMC have 
adopted the IAU standard values of d = 8.5 kpc and 
Vdrc = 220 km s _1 , as derived by Kerr & Lynden-Bell 
(1986; hereafter KLB86). These values assume the angu- 
lar velocity of circular rotation for the Sun to be 0o/d = 
25.9 km s _1 kpc _ . RB04 directly measured this valuqjfrom 
the proper motion of Sgr A* to be 0o/d = 29.45 ± 
.15 km s _1 kpc _1 , which using d = 8±.5 kpc (Reid 1993) 
implies Vdrc = 236±15 km s _1 . Based on the observed or- 
bits of individual stars around Sgr A*, Ghez et al. (2008) and 
GGE08 measured the distance of the Sun from the MW cen- 
ter to be d = 8.4±0.4 kpc, assuming that Sgr A* is at res10. 
These latest measurements are similar to the value inferred 
by KLB86, and bring the maximum Vdrc (within 1-a) up to 
265 km s _1 , which is ~ 20% higher than the 220 km s _1 on 
which all prior studies were based. Calculations based solely 
on KLB86 are therefore at the lower end (2-cr below) of 
the allowed values for M V i r , Vdrc, and r V i r - The new values 
are listed in Table [T] The inferred masses and radii assume 
the commonly accepted value of M V i r = 1 x 1O 12 M for 
Vdrc = 220 km s _1 . The values of vx, W, vz (correspond- 
ing to Galactic Coordinates X, Y, Z) were all calculated us- 
ing the standard method of vdM02, changing only the value 

of Vdrc- 

While all of the recalculated masses are consistent 
with various recent measurements, the circular velocity of 
251 km s _1 , corresponding to M„ ir = 1.485 x 1O 12 M , gives 
the minimum value of the virial mass that is consistent with 
the timing argument, which puts the total mass of the Lo- 
cal Group (LG) between 3.2 x lO 12 Af and 5.5 x 1O 12 M 
(Binney & Tremaine 1987, hereafter BT87; van der Marel & 
Guhathakurta 2008, hereafter vdMG08; Li & White 2008). 
The MW and M31 galaxies dominate the LG and are of 
comparable masses; a combined mass of 3-4 x 1O 12 M is 
on the lower end of the timing argument estimate. Li & 
White (2008), using the Millennium Simulation and data 
on Leo I, find the most likely mass for MW halo to be 
M M w = 2.34 x 1O 12 M with a lower limit of 0.8 x 1O 12 M , 
giving a more than adequate range to acommodate our in- 
ferred masses. Additional mass estimates of the MW halo 
are found in Table [3] and discussed further in 35] 

The P08(251 km s^ 1 ) entry (with this notation de- 
noting the proper motion from P08 and a circular veloc- 
ity of 251 km s" 1 ) in Table 2 for the LMC velocity of 
339 km s _1 is significantly lower (by more than 2-a) than 
the K06(220 km s _1 ) value of 378±18 km s _1 , although it 
is higher than the vdM02(220 km s _1 ) weighted average 



1 The supermassive black hole, Sgr A*, is expected to lie nearly 
motionless at the dynamical center of the MW since its mass far 
exceeds that of the surrounding stars. Hence, the apparent motion 
of SgrA* on the sky represents the reflex of the Sun's orbit around 
the Galactic center, including both the mean Galactic rotation 
and the small peculiar motion of the Sun relative to the local 
standard of rest (RB04). 

2 Since the surrounding stars are lighter by six orders of magni- 
tude than Sgr A*, and since there is currently no evidence for a 
second (intermediate mass) black hole, the assumption that Sgr 
A* is at rest at the dynamical center of the MW appear most 
natural. 
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Table 1. Modified Values for V c i rc , M v i r , and r v i r of the MW 









color 




(kpc) 


(km s — 1 ) 


(10 12 Af Q ) (kpc) 






7.5 


220 


1.000 258 


yellow 




8.0 


236 


1.234 277 


green 




8.5 


251 


1.485 295 


blue 




9.0 


265 


1.748 310 


purple 




Table 2. Calculated velocities for the LMC 




Author 




Vx , Vy , vz 


V 


color 




(km s- 1 ) 


(km s" 1 ) 


(km s _1 ) 




vdM02 


220 


-56, -220, 186 


293 




vdM02 


251 


-56, -189, 186 


271 


gray 


K06 


220 


-86, -268, 252 


377 




K06 


251 


-86, -237, 252 


356 


black 


P08 


220 


-83, -258, 238 


360 


yellow 


P08 


236 


-83, -243, 238 


350 


green 


P08 


251 


-83, -227, 238 


339 


blue 


P08 


265 


-83, -213, 238 


330 


purple 



of previous studies, 293±39 km s" 1 . Comparing the cor- 
responding LMC velocities for Vcirc = 251 km s _1 , the 
K06 and P08 values are much closer (in agreement with 
their initial transverse velocities being within one stan- 
dard deviation of each other). P08, K06, and vdM02 with 
Vdrc = 251 km s _1 are all compared in Figured] Figure [2] 
compares the P08 proper motion at various values of Vdrc- 
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Figure 1. Distance of the LMC from the center of the MW galaxy 
as a function of time (with zero as the present time) for orbital 
velocities published in previous studies (see Table 2), but assum- 
ing V c%T c = 251 km s" 1 for the MW. The blue line is from P08, 
the gray line is from vdM02, and the black line is from K06/B07. 
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Figure 2. Distance of the LMC from the center of the 
MW galaxy. The yellow trajectory (corresponding to V c irc = 
220 km s _1 ) is gravitationally unbound to the MW. The green 
trajectory (corresponding to Vdrc = 236 km s — 1 ) is also un- 
bound. The blue orbit (corresponding to V c irc = 251 km s _1 ) is 
bound to the MW, with a period of ~ 6.3 Gyr and an apogalac- 
ticon of 347 kpc. The value of r v i r at this Vdrc is 295 kpc 
(see Table [l] and equation O . The purple orbit, corresponding 
to Vdrc = 265 km s — x , is more tightly bound, with a period of 
~ 3.5 Gyr and an apogalacticon distance of 214 kpc. This dis- 
tance is well within the virial radius r v i r = 310 kpc at this value 

of Vcirc 



4 TIDAL EFFECT OF M31 

In our simplified analysis of the dynamics of M31, we ig- 
nore its transverse motion and consider only the perturba- 
tive influence of M31 as its radial distance changes relative to 
the MW (BT87). While the transverse component of M31's 
velocity might be non-negligible (Loeb et al. 2005), recent 
analysis suggests that it is lower than the radial component 
(vdMG08) and its inclusion would only have a weak effect on 
the results reported here. We also ignore any diffuse inter- 
galactic mass in between the MW and M31, although future 
modeling might take the related uncertainty into considera- 
tion (Cox & Loeb 2007). To account for the changing sepa- 
ration between the two galaxies, we used the standard radial 
dynamics model (BT87). 

We find that M31 has a small but non-negligible effect 
on the trajectory of the LMC. Based on the current position 
of the LMC (see Figures |3] and gj, the addition of M31 pulls 
the LMC away from the Galactic center and towards the LG 
center of mass. Figure [5] compares the distance of the LMC 
from the Galactic (MW) center with (blue line) and without 
(light blue line) the gravitational influence of M31, assuming 
Vdrc = 251 km s" 1 for the MW. Our figures keep the MW 
at the center of the coordinate system (corresponding to an 
accelerated frame of reference). Without M31, the LMC 
is more tightly bound to the MW; in the P08(251km s" 1 ) 
case, the orbital period of the LMC is 5 Gyr, about 1.5 Gyr 
shorter than if M31 is included. 
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Y (kpc) 



Figure 3. Projection of the path of the LMC onto the Y— Z 
plane (Galactic coordinates) of the LG over the past 8 Gyr. 
The yellow, green, blue, and purple trajectories correspond to 
V circ = 220, 236, 251, and 265 km s _1 , respectively. The red 
line traces the path of M31 (the large red dot) from and towards 
the MW. The Galactic (MW) center is located at the black dot 
and the LMC center is represented by the small blue dot. The 
approximate virial radii of the MW and M31 are traced out by 
the dashed lines. 

z (kpc) 




Figure 4. A close-up of the MW galactic center in Figure [3] 
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Time (Gyr) 

Figure 5. Distance of the LMC from the center of the MW as a 
function of time, with (darker line) and without (lighter line) the 
gravitational influence of M31, assuming V c i rc = 251 km s — 1 . 



5 DISCUSSION 

We have found that the ~ 14 ± 6% increase in the MW 
circular velocity, relative to the IAU standard of 220 km s _1 , 
allows the LMC to be gravitationally bound to the MW. 
Despite its relatively high proper motion (K06; P08), the 
orbit of the LMC remains confined within the virial radius 
of the MW. 

Since the MW and M31 galaxies account for most of the 
mass in the LG, the LG itself is not much more spatially ex- 



tended than the two are individually, having an estimated 
zero energy surface at a radius of ~ 0.9 Mpc (Karachent- 
sev et al. 2008). If the LMC had followed the path dic- 
tated by the P08(220 km s" 1 ) parameters in Table 2, it 
would have originated from a distance of > 1.5 Mpc away 
from the MW center (see Figure |3J). In comparison, if we 
admit the P08(251 km s _1 ) parameters, then the LMC 
originated roughly at the virial radius of the MW halo 
and well within the boundaries of the LG. The path sug- 
gested by P08(236 km s" 1 ) puts the LMC origin outside 
the LG but closer than the P08(220 km s _1 ) case. If the 
P08(265 km s _1 ) model is to be believed, the LMC is on its 
third pass by the MW center and also originated on the edge 
of the MW halo. The Galactocentric distance on which this 
result is based (d© = 9.0 kpc) is 1-a above the most recent 
estimate (8.4T0.5 kpc; see GGE08), so it is not excluded. 

We note that during the long time-scale between perigee 
passes (~ 6 Gyr, see Fig. [2j , the MW halo mass could evolve 
by tens of percent (Diemand, Kuhlen & Madau 2007). Fu- 
ture studies might use cosmological simulations to incorpo- 
rate the evolution of both the host (MW-like) galaxy and 
its most massive (LMC-like) satellite to get a statistical un- 
derstanding of their likely past interaction. Our most likely 
model has the LMC forming at roughly the MW virial ra- 
dius, as expected from a hierarchical formation of the MW 
halo (Springel et al. 2008, Figs. 11 & 12). Also, all our quoted 
uncertainty for the LMC orbit stem from the uncertainty in 
d (and therefore Vcirc)- Errors in the proper motion were 
not taken into account. The K06 measurements, for exam- 
ple, are slightly over 1-a off from the P08 measurements of 
the LMC proper motion. The difference between these two 
cases (as seen in Figure [TJ is drastic - in one case the LMC 
is clearly bound to the MW (P08) , and in the other it is not 
(K06). A 1-a shift in the other direction, however, would 
alter the path of the LMC into an even tighter orbit, similar 
to the P08(265 km s" 1 ) model. 

With the value of the distance between the Sun and 
the Galactic center at its IAU value (GGE08, KLB86), 
~8.5±.5 kpc, the total mass of the MW Galaxy increases 
by a factor of 1.23-1.75 relative to the values inferred for 
the lower distance of 8±.5 kpc (Reid 1993). B07 correctly 
rules out a larger mass for the MW but only considers the 
low mass (1 x lO 12 A/ ) and the high mass (2 x 1O 12 M ) 
models of Klypin et al. (2002) which bracket our preferred 
range. Xue et al. (2008) describe the uncertainties of the 
prior assumption that the galactic Vcirc is 220 km s _1 . Ta- 
ble [3] provides the corresponding range of masses for the MW 
halo, all calculated from observational data, fit either to a 
flat rotation curve or an NFW profile, such as the one we 
have adopted in this work. 

The ~ 50% increase in mass (from 1 x 10 12 M to 
1.485 x 10 12 M Q ) and the ~ 6% decrease in the velocity 
of the LMC (P08(220 km s" 1 ) to P08(251 km s" 1 )) have a 
comparable influence on making LMC orbit bound to the 
MW. The combined effect of these changes is not equivalent 
to a change in the MW mass, as considered by B07. The 
K06(251 kms -1 ) and P08(220 km s" 1 ) lines in Figures [T] 
and [2] correspond to roughly the same LMC velocity (as de- 
rived from the IAU standard with the P08 proper motion) 
but a different MW mass. The higher mass of the MW con- 
tributes about half of the overall difference between these 
cases. This parallels the comparison made in B07 between 
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Table 3. Milky Way Halo Mass in Recent Studies 



Author 


Mass 


Model 




(1O 12 M ) 





Wilkinson & Evans (1999) 
Sakamoto, Chiba & Beers (2003) 
Sakamoto, Chiba & Beers (w/o Leo I) 
Smith et al. (2007) 
Xue et al. (2008) 

Li & White (2008) 



1 9+ 3 ' 6 
^+0.5 
-1.0 

1.8j- M 



2.5: 



-0.7 
,+1.14 

-0.54 

0.79±0.15 - 



1.42" 1 



1.18 ±0.28 
2.43 



FRC 1 
FRC 
FRC 
NFW 
NFW 

N-Body 2 



X FRC denotes a Flat Rotation Curve, as in an isothermal sphere. 
These mass calculations have cutoffs of ~ 50 kpc. 2 e.g. the 
Millennium simulation. Here the limiting radius is the virial 
radius, similar to the NFW calculation. 

the K06(220 km s _1 ) and the K06(220 km s" 1 )+4cr num- 
bers in their Fiducial and High Mass Models, except that 
our findings do not require a 4 — a deviation from the mea- 
sured LMC velocity or a 2 — a deviation from the mass of 
the MW in order to make the LMC orbit bound to the MW. 

The previous suggestion (B07) of possibly disqualifying 
the LMC and SMC as MW satellites has undesirable impli- 
cations for the satellites in the MW halo. M31 has 18 satel- 
lites, 5 of which are gas rich, whereas the MW, if the LMC 
and SMC are no longer bound to it, has only 12 bound satel- 
lites, 2 of which are gas rich (Karachentsev 2005, vdMG08). 
M31's satellites range in mass from 0.58 to 500 xl0 8 M©, 
whereas the remaining MW satellites are in the range of 
0.1-1 xlO 8 M (Mateo 1998). Excluding the LMC and SMC 
as satellites of the MW appears unnatural, as there is no 
reason for the two comparably sized galaxies to have a large 
disparity in the abundance of massive satellites. Moreover, 
the chance of finding massive galaxies like the LMC and the 
SMC so close to the MW center requires a special coinci- 
dence if they are unbound to the MW, since they would 
have spent most of their orbital time far away from the MW 
in that case. Yet, no similar galaxies are known to exist 
in the much larger volume between M31 and the MW. Al- 
though these statistical arguments are not definitive, they 
support indirectly the higher updated values of V c i rc (RB04) 
and d (GGEO8) for the MW. 

No work on the Magellanic Clouds would be complete 
without a mention of the Magellanic Stream (MS). There 
is < 2% change in the trajectory over the length of the 
MS (100 degrees) between the P08(220km s" 1 ) and the 
P08(251km s" 1 ) models, both of which are well within the 
error margins in Figure 8 of B07. We reiterate the concerns 
presented in B07 that neither tidal stripping nor ram pres- 
sure can fully explain the orientation of the Stream. 
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